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Abstract 



We present an analysis of the neutrino oscillation solutions of the solar and at- 
mospheric neutrino problems in the framework of four-neutrino mixing where 
a sterile neutrino is added to the three standard ones and the mass spectra 
presents two separated doublets. Such scenarios allow for simultaneous tran- 
sitions of solar z^e, as well as atmospheric v^, into active and sterile neutrinos 
controlled by the additional mixing angles 1)23 and ??24, and they contain as 
limiting cases the pure solar fe -active and -sterile neutrino oscillations, 
and the pure atmospheric Vg and oscillations, respectively. We 

evaluate the allowed active-sterile admixture in both solar and atmospheric 
oscillations from the combined analysis. Our results show that, although the 
Super-Kamiokande data disfavour both the pure — > atmospheric chan- 
nel and the pure v^. — > Vg solar channel, the result from the combined analysis 
still favours close-to-pure active and sterile oscillations and disfavours oscil- 
lations into a near-maximal active-sterile admixture. 
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I. INTRODUCTION 



Super-Kamiokande high statistics data indicate that the observed deficit in the 
/i-hke atmospheric events is due to the neutrinos arriving at the detector at large zenith 
angles, strongly suggestive of the oscillation hypothesis. Similarly, their data on the 
zenith angle dependence and recoil energy spectrum of solar neutrinos in combination 



with the results from Homestake p, SAGE [|], and GALLEX+GNO |TO|,|ni experiments, 
have put on a firm observational basis the long-standing problem of solar neutrinos, strongly 
indicating the need for conversions. 

In addition to the solar and atmospheric neutrino results from underground experiments, 
there is also the indication for neutrino oscillations in the Ue channel by the LSND 



experiment [|T2|. All these results can be accommodated in a single oscillation framework 
only if there are at least three different scales of neutrino mass-squared differences. The 
simplest case in which this condition is satisfied requires the existence of a fourth light 
neutrino, which must be sterile {i.e. having interactions with standard model particles much 
weaker than the SM weak interaction) in order not to affect the invisible Z decay width, 
precisely measured at LEP [pr^-p!7[]. 

There are various analyses in the literature comparing the oscillation channels into ac- 
tive or sterile neutrinos for both solar PT|-P3| and atmospheric data p!H| , P!^ performed in 
the framework of two-neutrino oscillations (where oscillations occur into either an active or 
a sterile state). In particular the Super-Kamiokande collaboration has claimed the oscilla- 
tion into sterile neutrinos to be disfavoured both for solar and atmospheric neutrinos [0. 
However, when considered in the framework of four-neutrino mixing, oscillations into pure 
active or pure sterile states are only limiting cases of the most general possibility of oscilla- 



tions into an admixture of active and sterile neutrinos |17]. Such a more general framework 



of four-neutrino mixing has also been analysed in recent studies, which have put constraints 



on the active-sterile admixture either in the context of solar p4| or atmospheric [p5| -|27 
neutrino oscillations, respectively. 

In this paper we present a combined analysis of the neutrino oscillation solutions of 
both the solar and the atmospheric neutrino problem, in the framework of four-neutrino 
mixing. We include in our analysis the most recent solar neutrino rates of Homestake [§, 
SAGE [0, GALLEX and GNO as well as the recent 1258 day Super-Kamiokande 

data sample 0, including the recoil electron energy spectrum for both day and night periods. 
As for atmospheric neutrinos we include in our analysis all the contained event from the latest 
79.5 kton-yr Super-Kamiokande data set [Q], as well as the upward-going neutrino-induced 
muon fluxes from both Super-Kamiokande and the MACRO detector ||28||. The constraints 
arising from relevant reactor (mainly from Bugey [EO]) and accelerator (CDHSW PD[ and 



CCFR |3T1]) experiments are also imposed. 

The outline of the paper is the following. In Sec. |I| we summarize the main expressions 
for the neutrino oscillation formulas that we use in the analysis of solar and atmospheric 
neutrino data which take into account matter effects both in the propagation in the Sun and 
in the Earth. Sections |IT1| and |IV| contain the results for the analysis of the four-neutrino 
oscillation parameters for solar and atmospheric data respectively. The results for the full 
combined analysis are described in Sec. |V|, with particular emphasis on the active-sterile 
oscillation admixture which better allows for a simultaneous description of the solar and 
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atmospheric neutrino data. Our results show that, against what may be naively expected, 
although the Super-Kamiokande data disfavours both the pure atmospheric channel 

and the pure solar channel, the result from the combined analysis still favours any 

of these close-to-pure active and sterile oscillations and disfavours oscillations into a near- 
maximal active-sterile admixture. 



II. FOUR-NEUTRINO OSCILLATIONS 

There are six possible four-neutrino schemes that can accommodate the results from 
solar and atmospheric neutrino experiments as well as the LSND evidence. They can be 
divided in two classes: 3 + 1 and 2 + 2. In the 3 + 1 schemes there is a group of three 
neutrino masses separated from an isolated one by a gap of the order of 1 eV^ which is 
responsible for the short-baseline oscillations observed in the LSND experiment. In 2 + 2 
schemes there are two pairs of close masses separated by the LSND gap. 3 + 1 schemes 



are disfavoured by experimental data with respect to the 2 + 2 schemes |I5|JT6|] , but they 
are still marginally allowed Therefore we concentrate in the 2 + 2 schemes shown in 
Fig. |T[ For the phenomenology of neutrino oscillations these two schemes are equivalent up 
to the relabeling of the mass eigenstates (or equivalently of the mixing angles). Thus in 
what follows we will consider the scheme B where the mass spectrum presents the hierarchy: 

Ami = ^^ll < ^^atm = ^^Is < A'^LSND = ^^ll - ^"^42 - ^"^31 - ^"^32 (2-1) 

(we use the common notation Am|j = m| — m^). In this four-neutrino mixing scheme, the 
flavour neutrino fields i^a (<^ = e, s, fi, r) are related to the mass eigenstates fields z/^ by 

4 

i^a = '^Uaki^k {a = e,s, fi,T) , (2.2) 
fe=l 

where is a 4x4 unitary mixing matrix. Neglecting possible CP phases, the matrix U can 
be written as a product of six rotations, Uu, U13, Uu, U23, f/24 and t/34 where 

{Uij)ab = Sab + {Cij - 1) (SiaSib + SjaSjb) + Sij {5ia6jb - ^ja^ib) ■ (2.3) 

where Cij = cos 'dij and Sij = sin-i^jj. The order of the product of the matrices corresponds 
to a specific parametrization of the mixing matrix U. In order to study oscillations of the 
solar and atmospheric neutrinos which include the matter effects in the Sun and/or in the 
Earth, it is convenient to use the following parametrization |T^: 

U = U24U23UuUrsU34Uu. (2.4) 
This general form can be further simplified by taking into account the negative results from 



the reactor experiments, in particular the Bugey one [29| which implies that 



in the range of Am^^ relevant for the LSND experiment. This leads to an upper bound on 
the projection of the over the 3-4 states: 
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e4 



(2.6) 



so that the two angles '(?i3 and must be small and their contribution to solar and at- 
mospheric neutrino transitions is negligible. Therefore, in what follows we will set these 
two angles to zero, although one must keep in mind that at least one of them must be 
non- vanishing in order to account for the LSND observation. 

In the approximation i^ia = di4, = the U matrix takes the form: 



U 



( Ci2 
-S12C23C24 
■Sl2"S23 
V -512023^24 



Sl2 
C12C23C24 
~Ci2S23 
"~Ci2C23S24 





S23C24C34 - S24S34 

C23C34 

~"S23S24C34 — C24S34 



\ 

S23C24S34 + S24C34 

C23"534 
~ 523524 "534 + C24C34 / 



(2.7) 



In the rest of this section we will discuss the consequences of the mixing structure of Eq. 
on the relevant transition probabilities for solar and atmospheric neutrino oscillations. The 



transition probabilities that take into account matter effects have been derived in Ref. [17 



Some improvement concerning the calculation of the regeneration of solar v^s in the Earth 
was presented in Ref. p3| . Here we summarize the discussion on the flavour composition of 



the relevant states for both solar and atmospheric neutrino oscillations. 



A. Solar Neutrinos 

In the scheme here considered, solar neutrino oscillations are generated by the mass- 
squared difference between V2 and Vi. It is clear from Eq. ( |2.?| ) that the survival of solar t'e's 
mainly depends on the mixing angle "(912, whereas the mixing angles '(923 and ^24, determine 
the relative amount of transitions into sterile Vs oi active i/^ and Vr- Let us remind the 
reader that and Vr cannot be distinguished in solar neutrino experiments, since their 
matter potential and their interaction in the detectors are due only to neutral-current weak 
interactions and therefore they are equal. Thus, instead of z/^ and z/,-, one can consider the 
linear combinations 




cos d sin d 
— sin -d cos d 




(2.8) 



with 



cos-d = , = , sin'i9 = 24C23 _ ^ {l.'i)) 



^23^24 V -L C23C24 

The mixing of Va and with Pi and P2 is given by 



Ual = -Suyl- Ci3ci4 , Ua2 = Ci2\J I - cl^cl^ , Ubl = Ub2 = , (2.10) 

so that the solar neutrino oscillations occur only between the states 



Pe^^a with Pa = C23C2iJ^s + y i - clsCl^Pa , (2-11) 

with mixing angle 1912, and 
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4Aa = 1 - WalV - \Ua2V = WslV + \Us2V , (2.12) 

gives the size of the projection of the sterile neutrino onto the state in which the solar u^. 
oscillates to. 

Therefore, solar neutrino oscillations depend only on ^9i2 and the product C23C24. We 
distinguish the following limiting cases of Eq. (|2.11|) : 

• if C23C24 = then Usi = Us2 = 0, Uai = — S12, Ua2 = cu, corresponding to the limit of 
pure two-generation t'e Va transitions; 

• if C23C24 = 1 then Usi = S12, Us2 = C12 and Uai = Ua2 = and we have the limit of 
pure two-generation Pg transitions. 

Since the mixing of Ve with vi and V2 is described by a two-generation equation with 
mixing angle di2-, the mixing of Vg with Vi and V2 is given by the same formula multiplied 
by C23C24, and the mixing of Va with vi and z/2 is again described by the same equation 



times Y 1 — C23C24, it is clear that, in the general case of simultaneous z/g — and z/g Va 
oscillations, the corresponding probabilities are given by 

pSun+Earth ^ ^2^^2^ _ pSun+Earth) ^ (2.13) 

= (i - 44) (i - C-^i^f^"') ■ (2-14) 

where P^^^yf^^^'^ takes the standard two-neutrino oscillation form for Am^i and 'du, 



found by numerically solving the evolution equation in the Sun and the Earth matter with 
the modified matter potential 

A = Acc + cIscI^Anc- (2.15) 

with Acc = 2y/2GFENe and Anc = —V^GpENn. Here A^e and Nn are, respectively, the 
number densities of electrons and neutrons in the medium, E is the neutrino energy and 
Gp is the Fermi constant. These expressions satisfy the relation of probability conservation 

pSun+Earth i pSun+Earth i pSun+Earth i 

The analysis of the solar neutrino data in the four-neutrino mixing schemes is therefore 
equivalent to the two-neutrino analysis but taking into account that the parameter space 
is now three-dimensional (Amg^, tan^ '!?i2, C23C24). Although originally this derivation was 
performed in the framework of the 2 + 2 schemes |T7| , p^ , it is equally valid for the 3 + 1 



ones 



Let us finally comment on the range of variation of the mixing parameters. Choosing 
the convention that Amgi > 0, from Eqs. (|2.13|) and ( p. 14 ) it is clear that for solar neutrino 



oscillations the full physical space is covered by choosing the mixing angles in the ranges: 

0<1?12<^, 0<l923<^, 0<1?24<^. (2.16) 

Actually, since the dependence on -(923 and ^24, enters only through the combination C23C24, 
for what concerns the solar neutrino phenomenology one of these mixing angles could be 
fixed to zero. However, as we will see in the next section, they enter independently in the 
relevant probabilities for atmospheric neutrinos. 
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B. Atmospheric Neutrinos 



Atmospheric neutrino oscillations are generated by the mass-squared difference between 
1^3 and z^4. From Eq. ( p.7| ) one sees that the atmospheric v^s are not affected by the four- 
neutrino oscillations in the approximation di-^ = = and neglecting the effect of Am^i 
in the range of atmospheric neutrino energies. Conversely the survival probability of at- 
mospheric z/^'s mainly depends on the mixing angle 'du- The mixing angles ^923 and '!?24 
determine the relative amount of transitions into sterile or active More explicitly 
Eq. (|2.7|) implies that the atmospheric neutrino oscillations, i.e. oscillations with the mass 



difference Arnl^ and mixing angle '(?34, occur between the states 

1^/3 — > with 1^/3 = S23C2il^s + C23t'^ - S23S34l^r and = S24l^s + C2aI't ■ (2-17) 

We distinguish the following limiting cases of Eq. ( |2.17| ) : 



• if C23 = 1 then U^i = U^2 = 0. The atmospheric ly^ = lyp state oscillates into a state 

= C24Z^T + 524 i^s. This is the limit studied in Ref. p5| , |26| . We will denote this case 
as "restricted" . In particular: 

— in the case C23 = C24 = 1 we have Us3 = Ug^ = 0, corresponding to the limit of 
pure two-generation ly^ — > 1/^- transitions; 

— in the case C23 = 1 and C24 = we have Urs = Uri = 0, corresponding to the limit 
of pure two-generation z/^ Ug transitions; 

• if C23 = there are not atmospheric neutrino oscillations as the projection of i/^ over 
the relevant states cancels out {U^s = U^^ = 0). 

Thus the mixing angle '323 determines the size of the projection of the over the "atmo- 
spheric" neutrino oscillating states, 

4 = + K2\' = 1 - Ksl' - (2.18) 

As a consequence, one expects S23 to be small in order to explain the atmospheric neutrino 
deficit and, as we will see in Sec. |V|, this is the case. Furthermore, the negative results from 
the CDHS and CCFR experiments on searches for t'^-disappearance also constrain 
such projection to be 



'23 



\U^,i\^ + \U^2\^ <0.2 (2.19) 



at 90% CL for Am^gNo ^0.4 eV^. 

Approximate expressions for the relevant survival probability including matter effects 
in the Earth can be found in Ref. [|1^]. For what concerns the work presented here, the 



relevant probabilities can be calculated by numerically integrating the evolution equation in 
the Earth with the modified matter potential 

A = (4 - s1^cI,)Anc , (2.20) 



so that for pure atmospheric — > z/^ oscillations (533 = -§24 = 0) ^4 = 0, while for Ug 
(^23 = 0, S24 = 1) A = A]sfc, as expected. It is this modification of the Earth matter 
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potential what gives to the atmospheric neutrino data the capabihty to discriminate between 
the active and sterile oscillation solution. In particular, higher sensitivity to this potential 
effect is achievable for the higher energy part of the atmospheric neutrino flux which lead 
to the upward going muon data. The main effect of the presence of this potential is that for 
pure active oscillations the angular distribution of upgoing muons is expected to be steeper 
at larger arrival angles, while a flattening is expected due to the matter effects for sterile 
oscillations P(]| . 

From Eqs. ( |2.17|) and ( |2.2CI| ) we see that, unlike in the case of solar neutrinos, the angles 
't?23 and t?24 enter independently in the atmospheric oscillations. Thus the analysis of the 
atmospheric neutrino data in the four-neutrino mixing schemes is equivalent to the two- 
neutrino analysis but taking into account that the parameter space is now four-dimensional 

^23; ^24)- 

Concerning the range of variation of the mixing angles, the full parameter space relevant 
for atmospheric neutrino oscillation can be covered by choosing 'i?23 and '(?24 in the ranges 
indicated in Eq. ( |2.16| ), but, in general, the mixing angle 7934 must be allowed to take values 
in the interval 

- 2 ^ ^34 < - (2.21) 

(in the convention Aml^ > 0). In the limiting case '(?23 = 0, there is an additional symmetry 
in the relevant probabilities so that the full parameter space can be spanned by < < | , 
as expected since in this case we have the effective two-neutrino oscillation ^ v^. 



III. SOLAR NEUTRINO ANALYSIS 

We first describe the results of the analysis of the solar neutrino data in terms of Ve 
oscillations. We determine the allowed range of oscillation parameters using the total event 
rates of the Chlorine P], Gallium and Super-Kamiokande (corresponding to the 

1258 days data sample) experiments. For the Gallium experiments we have used the weighted 
average of the results from GALLEX+GNO and SAGE detectors. We have also included 
the Super-Kamiokande electron recoil energy spectrum measured separately during the day 
and night periods. This will be referred in the following as the day-night spectra data which 
contains 19 + 19 data bins. The analysis includes the latest standard solar model fluxes. 



BPOO model [^, with updated distributions for neutrino production points and solar matter 
density. 

For details on the statistical analysis applied to the different observables we refer to 



Ref. [pT] , p3| . As discussed in the previous section, the analysis of the solar neutrino data in 
these four-neutrino mixing schemes is equivalent to the two-neutrino analysis but taking 
into account that the parameter space is now three-dimensional (Am2]^, tan^ '!?i2, CI3CI4), so 
we have a total of 37 degrees of freedom (d.o.f.). 

We first present the results of the allowed regions in the three-parameter space for the 
global combination of solar observables. Notice that since the parameter space is three- 
dimensional the allowed regions for a given CL are defined as the set of points satisfying the 
condition 

XL(A'^L, ^12, 44) - XsoLmin < Ax'(CL, 3 d.o.f.) (3.1) 
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where, for instance, Ax^(CL, 3 d.o.f.) = 6.25, 7.81, and 11.34 for CL = 90%, 95% and 99% 
respectively, and Xsoi min is the global minimum in the three-dimensional space. 

In Fig. 1^ we plot the sections of such a volume in the plane (Amgi, tan^('i?i2)) for different 
values of the active-sterile admixture |f/sip + |t^s2p = c\.^c\^- The global minimum used in 
the construction of the regions lies in the LMA region and for pure i/g^active oscillations 

Am^i = 3.65 X 10'^ eV^ 
tan^ d2i = 0.37 
Xsoi,min = 37.1 (3.2) 

which for 37 d.o.f. (3 rates + 38 spectrum points — 1 spectrum free normalization — 3 
parameters) corresponds to a goodness of the fit (GOF) of 46%. 

As seen in Fig. ^ the SMA region is always a valid solution for any value of C23C24 at 99% 
CL. This is expected as in the two-neutrino oscillation picture this solution holds both for 
pure ^active and pure t'e^sterile oscillations. Notice, however, that the statistical analysis 
is different: in the two-neutrino picture the pure i/g^active and t'e^sterile cases are analysed 
separately, whereas in the four-neutrino picture they are taken into account simultaneously 
in a consistent scheme. On the other hand, both the LMA and LOW solutions disappear 
for larger values of the active-sterile admixture. 

Let us comment here that, unlike for the case of atmospheric neutrinos where the dis- 
crimination between active and sterile oscillations arises from the difference in the matter 
potentials, for solar neutrinos the main source of difference is due to the lack of neutral 
current contribution to the water Cerenkov experiments for the sterile case. Unlike ac- 
tive neutrinos which lead to events in the Super-Kamiokande detector by interacting via 
NC with the electrons, sterile neutrinos do not contribute to the Super-Kamiokande event 
rates. Therefore a larger survival probability for ®B neutrinos is needed to accommodate the 
measured rate. As a consequence a larger contribution from ®B neutrinos to the Chlorine 
and Gallium experiments is expected, so that the small measured rate in Chlorine can only 
be accommodated if no ''Be neutrinos are present in the flux. This is only possible in the 
SMA solution region, since in the LMA and LOW regions the suppression of ^Be neutrinos is 
not enough. Notice also that the SMA region for oscillations into sterile neutrinos is slightly 
shifted downwards as compared to the active case. This is due to the small modification 
of the neutrino survival probability induced by the different matter potentials. The matter 
potential for sterile neutrinos is smaller than for active neutrinos due to the negative NC 
contribution proportional to the neutron abundance. For this reason the resonant condition 
for sterile neutrinos is achieved at lower Am^. 

In Fig. ^we plot the difference Ax^oi between the local minimum of Xsoi GSich solution 
and the global minimum in the three-dimensional space, as a function of the active-sterile 
admixture \Usi\'^ + |?7s2p = c^^c^^^. From the figure we find that: 

• solar neutrino data favours pure Ue Va oscillations but sizable active-sterile admix- 
tures are allowed; 

• the dependence of Ax^^j on the active-sterile admixture is very gentle in the SMA 
region while it is much more pronounced for the LMA solution; 
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• the three dimensional regions disappear at 90% (99%) CL for the following value of 

„2 2 
•-23 '-24 

Wsil"^ + |f/.2p = 4^4^ < 0.08 (1.) for SMA, 

+ |f/s2p = c^3C24 < 0.43 (0.71) forLMA, (3.3) 
\Usi\'^ + |f/s2r = C23C24 < 0.3 (0.77) for LOW-QVO; 



• at 95% CL the three-dimensional LMA region is allowed for maximal active-sterile 
mixing C23C24 = 0.5 while at 99% CL all solutions are possible for maximal admixture. 

There is a subtlety on the meaning of these limiting values of the mixing C23C24. The values 
quoted in Eq. ( |3.3| ) correspond to the mixings at which the corresponding three-dimensional 
region disappears at a given CL. Strictly speaking this is not the same as the allowed range of 
C23C24 at a given CL when the parameters (tan^ 'di2, Am^i) are left free to vary, either in the 
full parameter space or within a given solution. Such limits would be obtained by imposing 

the condition = Xsol,mm 12(C23C24) -XsoLmin < ^^^(CL, 1 d.O.f.) where xL,mml2i4344) 

is minimized in the parameter space 1-2 (or in a given region of such parameter space) and 
Xsoimin is minimized in the same region of the space 1-2 and in the mixing C23C24. 

In what follows we label as "constrained" the results of analyses when the parameters 
Am^^ and '^12 are varied within a given solution while "unconstrained" refers to the case 
where we allow the variation of Arn^i and in the full plane. For the unconstrained 
fit the corresponding curve in Fig. § is simply the lower envolvent of the curves from the 
constrained analysis, and it follows the LMA curve for < C23C24 < 0.46 and the SMA one 
for 0.46 < 4.4, < 1. 



IV. ATMOSPHERIC NEUTRINO ANALYSIS 



In our statistical analysis of the atmospheric neutrino events we use all the samples of 
Super-Kamiokande data: e-like and /i-like data samples of sub- and multi-GeV P], each 
given as a 5-bin zenith-angle distribution[|, and upgoing muon data including the stopping 
(5 bins in zenith-angle) and through-going (10 angular bins) muon fluxes. We have also 
included the latest MACRO [^] upgoing muon samples, with 10 angular bins which is also 
sensitive to the active-sterile admixture. So we have a total 45 independent inputs. 

For details on the statistical analysis applied to the different observables we refer to 



Ref. 1 18,35]. As discussed in the previous section the analysis of the atmospheric neutrino 



data in these four-neutrino mixing schemes is equivalent to the two-neutrino analysis but 
taking into account that the parameter space is now four-dimensional (Amfg, 1934, C23, C24). 
Alternatively we present the results in the equivalent parameter space (Amig, '!934, + 

\Uf.2\MUsl\^ + \Us2\^). 

We first present the results of the allowed regions in the four-parameter space for the 
global combination of atmospheric observables. Notice that since the parameter space is 



*Note that for convenience and maximal statistical significance we prefer to bin the Super- 
Kamiokande contained event data in 5, instead of 10 bins. 
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4-dimensional the allowed regions for a given CL are defined as the set of points satisfying 
the condition 



Xltmi^ml,, ^934, 4, 4) - Xatm.min < ^x\CL, 4 d.O.f.) (4.1) 

where Ax^(CL, 4 d.o.f.) = 7.78, 9.49, and 13.3 for CL = 90%, 95% and 99% respectively, 
and Xatm min the global minimum in the four-dimensional space. 

In Fig. ^we plot the sections of such a volume in the plane (Am^g, sin^('!934)) for different 
values of the mixings and 'i?24, which we parametrize by values of the projections + 
|[/^2p and + |t^s2p- As discussed in Sec. y for arbitrary values of C23 and C24 the full 

parameter space is covered for '!?34 in the range given in Eq. (|2.21| ). In Fig. ^ we display 
this full parameter space by showing the region in (Am43, sgn('!?34) sin^(i934)) with "positive" 
sin^(t?34) for < < | and "negative" sin^i^^si) for — | < -^34^ < 0. As seen from the figure 
the parameter space is composed of two separated regions, each of them around maximal 
mixing 'd34 = ±f • 

The global minimum used in the construction of the regions lies almost in the pure 
atmospheric ly^-i^r oscillations. At the best fit point: 

|f/si|'+ 1^7,21^ = 44 = 0.97 
\U^i\'' + \U^2\^ = sl, = OM 

Am^3 = 2.4 X 10"^ eV^ 
^^34 = 39° 

Xatm,min = 29.0 (4.2) 

which for 45 — 4 = 41 d.o.f. corresponds to a GOF of 92%. A careful reader may find this 
GOF surprisingly high. Let us comment here that there has been a reduction of the 
with respect to previous analysis [|T^,§^ due to a better agreement of the 79.5 SK electron 
distribution with the non-oscillation expectation (same effect is found in Ref. @]). Since 



the scenarios here discussed do not affect the atmospheric i^e predictions, this implies an 
overall improvement on the GOF of any these four-neutrino oscillation schemes. 

From Fig. ^ we see that the region becomes considerably smaller for increasing values 
of the mixing angle "^23, which determines the size of the projection of the t'^ over the 
"atmospheric" neutrino oscillating states, and for increasing values of the mixing angle 
1^24, which determines the active-sterile admixture in which the "almost-//^" oscillates to. 
Therefore from the analysis of the atmospheric neutrino data we obtain an upper bound on 
both mixings which, in particular, implies a lower bound on the combination C23C24 — 
|C/s2p limited from above by the solar neutrino data. To quantify these bounds we display 
in Fig. ||.a the allowed region in the parameter space (s^g = + |f^/x2p5 44 — l^siP + 

|[/s2p)- Following the discussion below Eq. ( p.3| ) we obtain the allowed two-dimensional 
region by imposing the condition 

AXatm = Xatm,min 34(^23, ^24) " xLn < Ax'(CL, 2 d.O.f.) (4.3) 

where Xatm mm 34(''^23, '*924) is minimized with respect to Amfg and {^34^ and Xatmmm is the 
global minimum in the four-dimensional parameter space. The diagonal cut in Fig. |^.a 
gives the unitarity bound 
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If/^ll' + K2\^ + If^^ll' + \Us2\^ = If^rSr + |f/r4r < 1 (4.4) 

where we have used that in Eq. ( p.7|) t/gs = f^e4 = 0. As seen from Fig. ||.a, the analysis 
of the atmospheric data imposes a severe bound on the z/^ projection. For instance, the 
two-dimensional region extends only till 

4 = If^MiT + I^m2|' ^ 0.12 (0.16) (4.5) 

at 90% (99%) CL, which, as previously mentioned (see Eq. (|2.19|) ), is perfectly consis- 
tent with the bounds from the z/^-disappearance accelerator experiments CDHDW and 
CCFR |31| given in Eq. (pl9|) . 



In order to determine the impact of the small but possible deviation from zero of S23 we 
also study the "restricted" case of Ref. [ p5| , p^ in which we impose s^^ = lU^il"^ + |f/^2p = 0. 
In this case the parameter space for atmospheric neutrino oscillation is three-dimensional 
and the allowed regions for a given CL are defined as the set of points satisfying the condition 

xl^{Aml,^U, 4) - Xatm,min < ^X^CL, 3 d.o.f.) . (4.6) 

In Fig. ^ we plot the sections of this volume in the plane {Aml^, sin^('(?34)) for different values 
of C24. As discussed in Sec. |I| for the "restricted" case the full parameter space is covered 
with < ^34 < |. This is graphically displayed in Fig. |, where the symmetry 'd^^ — > — 1934 
is clearly evident. 

The constrains on the active-sterile admixture in the two cases are displayed in Fig. |^.b, 
where we plot the dependence of Axltm on the active-sterile admixture C23C24 = 
|f^s2p5 both in the general case {i.e. when Xatm minimized with respect to all other pa- 
rameters in the problem), and in the "restricted" case (for which the projection has been 
fixed \U^i\^ + |^7^2p = 0). The horizontal lines are the values of Ax^(CL, 1 d.o.f.) for CL = 
90% and 99%. From this figure we obtain the 90% (99%) CL lower bounds on C23C24 from 
the analysis of the atmospheric neutrino data: 

C23C24 = \Usi\^ + \Us2\^ > 0.64 (0.52) (4.7) 
C24 = + |f/.2r > 0.83 (0.74) for the restricted case. (4.8) 

As seen in the figure Axltm has a monotonically growing behaviour as |^7sip + |t/s2p decreases 
for |^7sip + |f/s2p ^ 0.4-0.6, while the dependence fiattens below that value with the presence 
of a shallow secondary minimum around |f/sip + |^s2p ~ 0.25-0.3. We have traced the 
presence of this secondary minimum to the combined effect of the behaviour of for the 
different sets of upgoing muon data. For all sets of upgoing muon data the best fit is very close 
to the global minimum in Eq. ( |4.2| ). MACRO data at the moment leads to a monotonically 
growing x^ as the sterile admixture grows while Super-Kamiokande upgoing data on both 
stopping and through-going muons presents a fiattening of their for [f/^ip + |?7s2p < 0.5- 
0.6. This lost of sensitivity of Super-Kamiokande at larger sterile admixtures arises from the 
fact that their latest upgoing muon data presents a fiattening of the muon suppression in the 
angular bins for arrival directions around —0.7 < cos 6 < —0.4 which is compatible with large 



sterile admixtures. The same effect has been found in Ref. [26|. When combining the x foi" 



Macro and Super-Kamiokande data samples and taking into account the strong correlation 
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between their respective theoretical errors, we find the appearance of this secondary shallow 
minimum. 

From the figure we see that oscillations into a pure sterile state are excluded at 99.97% 
CL as compared to pure active oscillations. One must notice, however, that taken as an 
independent scenario, Vg oscillations give a x^j„ = 46 (for 45 — 2 d.o.f.) which implies 

a GOF of 35% and therefore cannot be ruled out in the bases of its absolute quality. 

In summary, we see that the analysis of the atmospheric data implies that: 

• one of the two-neutrino states in the 3-4 oscillating pair must be a close-to-pure fj^; 

• this close-to-pure oscillates into a state which is preferably composed by al- 
though a 48% admixture of sterile in this state is still allowed at 99% CL; 

• when imposing the "restricted" condition that one of the states is a pure v^^ the bound 
on the sterile admixture is tightened approximately by a factor 2 (from 48% to 26% 
at 99% CL). 

V. COMBINED ANALYSIS: CONCLUSIONS 

From the previous sections we have learned that the analysis of the solar data favours 
the scenario in which the solar oscillations in the plane 1-2 are Vf. oscillations into an active 
neutrino, and from that analysis we found an upper limit on the projection of the Vg on the 
1-2 states. On the other hand, the atmospheric neutrino analysis prefers that the oscillations 
of the 3-4 states occur between a close-to-pure and an active (z/^) neutrino, so giving an 
upper bound on the projection of the over the 3-4 states or equivalently a lower bound 
on its projection over the 1-2 states. The open question is then what is the best scenario 
for the active-sterile admixture once these two bounds are put together. 

To address this question in full generality we have studied the behaviour of the global 
Xltm+soi function defined as 

Xatm+sol(^"^21> "^12, Am^g, ??34, 'i?23, '!^24) = 

Xsol(A"^21, "^^23, C23C^4) + Xatm(A^43' ^^34, C23, 0^4) (5.1) 

with the mixing angles ^?23 and 1^24 determining the active-sterile admixture for the following 
scenarios: 

• ATM+SOLunc^ in this case we minimize Xatm+soi with respect to Am^i, 't?i2, Am^g 
and ??34, allowing Am|j^ and "Qyi to change from LMA to SMA as the active-sterile 
admixture controlled by -diz and ■j?24 change. The projection — \Upx^ + \ Utx2^ 
is free to vary in the range allowed by the analysis; 

• ATMR+SOLunc: in this case we fix = + \Uija^ = 0, so \Usi\^ + |t/s2p = C24. 
We minimize Xatm+soi with respect to Am|i, '!?i2, Amfg and ■i?34, allowing Am|i and 
■j?i2 to change from LMA to SMA as the mixing angle ■j?24 changes; 

• ATM-I-SOLlma(sma): in this case we minimize Xatm+soi with respect to Am^^, ■i?i2, 
Amig and 't?34 with Arn^i and '(?i2 constrained to be in the LMA (SMA) region. The 
i/fj^ projection s^g = l^^^uiP + \ U^2\'^ is free to vary in the range allowed by the analysis; 
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• ATMr+SOLlma(sma): in this case we fix 5^3 = + \U^2\'^ = 0, so |f/sip + |t/s2p = 

. We minimize Xatm+soi with respect to Aml^, ^12, Aml^ and 'dsi, with Amgi and 
'di2 constrained to be in the LMA (SMA) region. 

In Table | we give the values Xmin.atm+soi of the best fit point and its associated GOF in 
the six-dimensional (five-dimensional for the "restricted" case) space in each of these cases 
together with the best fit value of the mixings |t/^ip + |f//i2p = S23 |f^siP + |f^s2p = C23C24. 
In Fig. 1^ we plot the shift in from the corresponding global minimum in each of the 
scenarios as the mixing |f/siP + |f^s2p = C23C24 varies. From the table and the figure we find 
the following behaviours: 

• there are two favourite configurations which we will denote as near-pure-sterile 
solar neutrino oscillations plus near-pure-active atmospheric neutrino oscillations 
(NPSS+NPAA), and close-to-active solar neutrino oscillations plus close-to-sterile 
atmospheric neutrino oscillations (CAS+CSA), respectively. NPSS+NPAA oscilla- 
tions are characterized by the minima 

+ I^s2p ~ 0.91-0.97 , (5.2) 



where the exact position of the minimum depends on the allowed admixture |f/^i| + 
|C/^2p- CAS+CSA oscillations are characterized by the minima 

|f/sir + |^:/s2r -0.18-0.2 . (5.3) 



In both the cases the minima are not very deep in (this is particularly the case for 
CAS+CSA); 

• in none of the cases maximal active-sterile (MAS) admixture |t/sip + |f/s2p = C23C24 = 
0.5 is favoured. 

The scenarios with the solar solution either within the SMA region (Fig. 0.c) or left uncon- 
strained (Fig. I^.a) prefer the NPSS+NPAA scenario. They give a better fit to the combined 
analysis with a GOF of 66%-68%. In these cases the dependence of Xatm+soi with the 
active-sterile admixture is more pronounced and it is dominated by the dependence of the 
atmospheric neutrino analysis, although there is a secondary minimum at the CAS+CSA 
configuration which for the ATM+SOLunc case is acceptable at 92% CL. In all these scenar- 
ios MAS admixture is ruled out at more than 99% CL. Qualitatively this behaviour arises 
from the fact that in these scenarios increasing the sterile admixture in the atmospheric pair 
is strongly disfavoured while it is still acceptable for the solar pair since the SMA solution 
is valid for pure —>■ solar oscillations. 

On the other hand, the scenarios with the solar solution within the LMA region (Fig. |^.b) 
prefer the CAS+CSA scenario, although the dependence of Xatm+soi with the active-sterile 
admixture is not very strong and it presents a secondary minimum near the complementary 
situation of close-to-sterile solar plus close-to-active atmospheric (CSS+CAA) scenario, 
which is acceptable at 90% CL. They give the worst fit to the combined analysis but they 
still present an acceptable GOF: 47.5% (56%) for the "restricted" (unrestricted) case. In 
these scenarios MAS admixture is ruled out at 99% CL for the "restricted" case and it is 
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acceptable at about 95% CL for the unrestricted one. This behaviour arises from the fact 
that in these scenarios increasing the sterile admixture both in the atmospheric pair and 
in the solar case is strongly disfavoured. It is precisely in this case that one would expect 
the combined solution to he near MAS admixture. However, the results show that it is the 
solar dependence (which dominates in the combination) what leads to the slightly better 
description for the CAS+CSA scenario with a secondary minimum at the complementary 
CSS+CAA scenario. 

Summarizing, in this paper we have performed a combined analysis of the neutrino 
oscillation solutions of the solar and atmospheric neutrino problem in the framework of 
four-neutrino mixing. Such scenarios allow for simultaneous transitions of solar as well 
as atmospheric into active and sterile neutrinos. Our results show that although the 
solar and atmospheric data analyzed independently disfavour the pure z/^ atmospheric 
channel and the pure v^. — > Vg solar channel, the result from the combined analysis still 
favours close-to-pure active and sterile oscillations either in the solar or in the atmospheric 
sector (depending on the selected solar solution) and disfavours oscillations into a near- 
maximal active-sterile admixture. 
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solar ^^^= atm 

LSND LSND 

FIG. 1. 2 + 2 mass schemes in four-neutrino mixing favoured from the combination of solar, 
atmospheric and LSND results 
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TABLE I. Xmin atm+sol ^est fit point and its associated GOF in the six-dimensional 

(five-dimensional for the "restricted" case) space in each of the cases discussed in the text together 
with the best fit value of the mixings = + If^^sP ^'^d C23C24 = \Usi\^ + |f^s2p- 
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FIG. 2. Results of the global analysis of solar neutrino data for the allowed regions in Am^x and 
tan^ ^12 for the four-neutrino oscillations. The different panels represent sections at a given value 
of the active-sterile admixture |?7<iiP + |f^s2p = C23C24 of the three-dimensional allowed regions at 
90%, 95% and 99% CL. The best-fit point in the three parameter space is plotted as a star. 
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FIG. 3. Difference Ax^ between the local minimum of for each solution of the global analysis 
of solar neutrino data and the global minimum in the plane (tan^ t?i2, Amg]^) as a function of the 
active-sterile admixture \Usi\'^ + |C^s2p = c\nC^^. 
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FIG. 4. Results of the analysis of atmospheric neutrino data for the allowed regions in Aml^ 
and 'i?34 for the four-neutrino oscillations. The different panels represent sections at given values 
of the projection + |?7^2p = S23 active-sterile admixture |C/siP + |?7s2p = ^2^(^24 

of the four-dimensional allowed regions at 90%, 95% and 99% CL. The best-fit point in the 
four-parameter space is plotted as a star. The last panel corresponds to the case in which has 
also been minimized with respect to the S23 and C23C24. 
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FIG. 5. Results of the analysis of atmospheric neutrino data for the allowed regions in 
Am43 ^^'^ ^34 four-neutrino oscillations in the "restricted" case of projection 

+ |f^A'2p = ^23 — ^- The different panels represent sections at given values of the ac- 
tive-sterile admixture |C/sip + |C^s2p = C24 of the three-dimensional allowed regions at 90%, 95% 
and 99% CL. The best-fit point in the three-parameter space is plotted as a star. 
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FIG. 6. (a) 90%, 95% and 99% CL allowed regions for the projection + |C/^2p = sis 

and the active-sterile admixture + |f^s2p = C23C24 from the analysis of the atmospheric data 

(see text for details), (b) A^^ as a function of the active-sterile admixture \Usi\^ + |f^s2p = C23C24. 
In the lower curve has been minimized with respect to 7934, Am43 pi'ojection 
|f^/iiP+|f^/i2p = sis- The upper curve corresponds to the "restricted" case | P + 1 C//i2 P = S23 — Oj 
for which has been minimized only with respect to ??34 and A?n,43. 
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FIG. 7. Ax^ as a function of the active-sterile admixture |f/sip + |f^s2|^ = C23C24. In 
the ATM curves have been minimized with respect to 1934 , Am43 ^'^'^ projection 
+ |f^/x2p = -sis' well as with respect to the 1-2 parameters t?i2 and Am2i. The ATMr 
curves correspond to the "restricted" case |C/^iP + |C^/i2p = ^23 — 0' ^o has been minimized 
only with respect to ^34, Am^g and the 1-2 parameters. In (a) the 1-2 parameters 'di2 and Am^j^ 
have been minimized in the full parameter space while in (b) [(c)] they have been restricted to lie 
in the LMA [SMA] solution region. 
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